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^ ! Abstract 

| On the cosmological length scale, recent measurements by WMAP have validated ACDM to a 



precision not see before in cosmology. Such is not the case on galactic length scales, however, where 
the 'cuspy-core' problem has demonstrated our lack of understanding of structure formation. Here, 
we propose a solution to the 'cuspy-core' problem based on the observation that with the discovery 



of Dark Energy, Aqe, there is now a universal length scale, X^e = c/i^-DEG) 1 ^ 2 , associated with 

the universe. This length scale allows for an extension of the geodesic equations of motion that 

affects only the motion of massive test particles; the motion of massless test particles are not 

' affected, and such phenomenon as gravitational lensing remain unchanged. An evolution equation 

^J" '_ 

for the density profile is derived, and an effective free energy density functional for it is constructed. 
We conjecture that the pseudoisothermal profile is preferred over the cusp-like profile because it 
has a lower effective free energy. A cosmological check of the theory is made using the observed 
rotational velocities and core sizes of 1393 spiral galaxies. We calculate as to be 0.68±o.nj this 



is within experimental error of the WMAP value 0.761±S"2S. We then calculate -R200 = 270±i3o 



kpc, which is in agreement with observations. We estimate the fractional density of matter that 
cannot be determined through gravity to be 0.197±o.oi7; this is nearly equal to the WMAP value 
for the fractional density of nonbaryonic matter 0.1961q'q26- The fractional density of matter that 
can be determined through gravity is then calculated to be 0.041^q'q3^; this is nearly equal to 
^ = 0.0416^|. 
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I. INTRODUCTION 

The recent discovery of Dark Energy (see and references therein) has broadened 

our knowledge of the universe, and has demonstrated once again that it can hold surprises 
for us. The discovery has, most assuredly, also brought into sharp relief the degree of our 
understanding of it. Only a small fraction of the mass-energy density of the universe is made 
up of matter that we have characterized; the rest consists of Dark Matter and Dark Energy, 
and the precise properties of either one is not known. They are nevertheless needed to explain 
what is seen on an extremely wide range of length scales. On the galactic (~ 100 kpc parsec), 
galactic (~ 10 Mpc), and supercluster (~ 100 Mpc) scales, Dark Matter has been used to 
explain phenomena ranging from the formation of galaxies and their rotation curves, to the 
dynamics galaxies and the formation galactic clusters and superclusters. On the cosmological 
length scale, Dark Matter and Dark Energy is needed to explain the observed evolution of 
the universe from the Big Bang to the present, and will determine its fate in the future. 

Observations thus tell us that over a vast range of length scales the dynamics and evolution 
of the observed universe is determined not by normal matter, but by Dark Matter and now 
Dark Energy (see [s| for a quantum-cosmology approach that does not require Dark Matter 
or Dark Energy). Yet while the need and invocation of Dark Matter is ubiquitous on a 
wide range of length scales, our understanding of how matter determines dynamics at the 
galactic length scale is lacking. Recent measurements by WMAP |4J have validated the 
ACDM model of cosmology to an precision not seen before in cosmology. The situation on 
the galactic scale is not nearly as settled, however. Here, the cuspy-core problem for the 
density profile of matter in galaxies , and 0, 0] for reviews) has demonstrated our 

lack of understanding of the formation of galactic rotation curves. 

Current understanding of the structure formation is based on the work of Peebles 10], 
where the seeds of galaxies are due to local fluctuations in the density of matter that grow as 
the universe expands. Analytical solutions of this model have been done [ll, 12, 13, 14] for 



a number of special cases, and have resulted in density profiles that are sensitive to initial 
conditions, and have a power-law dependence whose exponents vary over a range of values. 
More recently, numerical simulations 

Asa 

of galaxy formation have been done, and have 
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consistently resulted in density profiles with a cusp-like structure 

Pi m 

PSimulation - fr/Rtf^ + (r/Rs)^'^ 1 ' 

instead of the expected pseudoisothermal density profile. Here, pi is a density parameter, 
R s is related to the radius of the galactic core, and the exponents (a, /3, 7) take a range of 
values from (1.5,2, 1.5) for Moore et. al. [?| to (1.0,3.0, 1.0) for Navarro, Frenk and Wilson 
(NFW) Q| (see Q] for review). This lead Moore to state in [?| that cold dark matter fails 
to reproduce the galactic rotation curves for dark-matter-dominated galaxies, one of the 
key reasons that dark matter was proposed in the first place. Soon afterwards, de Blok 



and coworkers 
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171 | explicitly demonstrated that the NFW density profile does not 



fit the density profile observed for Low Surface Brightness (LSB) galaxies (see also [18j for 
a recent analysis of cusp structure). Rather, the traditional pseudoisothermal profile, with 
(a,/3, 7) = (2,2,0), is the better fit. This demonstration is especially compelling as it is 
believed that dark matter dominates dynamics in LSB galaxies. 

There have been a number of attempts to solving the cuspy-core problem within ACDM 



19 . 20 . 21 , 22 1, and they have had varying degrees of success (see (8[ for a review ) . This 
problem does not exist in Milgrom's Modified Newtonian Dynamics (MOND) 23|, |24], [25] 
(see [26(] for a review) — a theory where Dark Matter is not needed — but there are a number 
of theoretical and observational problems that MOND must overcome (see 23] for arguments 
in support of MOND, however). 

Our approach to solving the cuspy-core problem, and to structure formation in general, 
is much more drastic; therefore, its reach is correspondingly broader. It is based on the 
observation that with the discovery of Dark Energy, Ade, we now have a universal length 
scale, Xde = c/^AdeG) 1 ^, on hand jz^J]. The geodesic equations of motion (GEOM) — and 
thus the geodesic action — is no longer unique, and extensions of it through the introduction 
of functions of Rc 2 / KoeG can be made. While there have been attempts at proving that the 
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31], such proofs 



GEOM are the unique consequence of the Einstein field equations 
assume that the background metric remains fixed under the passage of the test particle. As 
our extension of the depends explicitly on the energy-momentum tensor of matter — which 
includes the motion of the test particle — these proofs do not preclude our extension of the 
GEOM. 



In form, our extension of the GEOM preserves the equivalence principal, and through the 
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choice of the function of Rc 2 /A^eG, we can insure that their effects are not measurable on 



Vsl) * 10 30 

g/cm? (from [4]) is so much smaller than any density of matter either presently achievable 



terrestrial scales. Physically, that this choice is possible is because Ade = (7.21^11) x 10" 



experimentally, or present in regions of space accessible to experiment. Correspondingly, 
Xde = 14020g ( j ) Q Mpc is much longer than the scale of any experiment that has been used to 
test general relativity. In fact, the issue is to make the theory relevant at galactic scales of 
a few kpc, and by doing so we arrive at an estimate for the exponent a\. This exponent is 
the only parameter in the theory, and it determines the power-law behavior of our extension 
of the GEOM. We also find that while affecting the motion of massive test particles, our 
extension does not affect the motion of massless test particles; photons still travel along null 
geodesies, and gravitational lensing and the deflection of light are left unchanged. 

Applied to galaxy formation, the extended GEOM reduces to a nonlinear evolution equa- 
tion for the density profile of a model galaxy in the nonrelativistic, linear gravity limit. This 
evolution equation minimizes a functional of the density, which is interpreted as an effec- 
tive free energy functional for the system. Solutions to this equation is found using various 
velocity curves for galaxies as driving terms, and these solutions are then used to calculate 
the free energy associated with various profiles. We conjecture that like Landau-Ginzberg 
phenomenological theories in condensed matter physics, the system prefers to be in a state 
that minimizes this free energy. Showing that the pseudoisothermal profile is preferred over 
cusp-like profiles reduces to showing that it has the lower free energy. 

In our model of a galaxy, the Hubble length scale Xh = c/H (where H = IiHq, is the 
Hubble constant, h = 0.732j£{jjg and H = 100 km/s/Mpc fl) naturally appears, even 
though a cosmological model is not mentioned either in its construction or in its analysis. 
What happens at galactic length scales are naturally tied to what happens at cosmological 
length scales with our approach; the combination of the Dark Energy length scale and the 
nonlinear aspects of the extended GEOM link the two. This linkage allows us to extrapolate 
from the statistical properties of the observed universe the properties of a representative 
galaxy. These properties are then used to provide a cosmological check of the theory. 

As with the Peebles model, the total density of matter for our model galaxy can be 
written as a sum of a background density Pasym P ( x ) and a linear perturbation p} J (x). As 
usual, Pasymp(x) does not contribute to the motion of stars within the galaxy, while Pj/(x) 
does; in the absence of forces other than gravity, observations of the dynamics of the stars 
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0.196tg;g| measured by 



will determine only p} 7 (x). But unlike the Peebles model, p a sym P is not a constant. With it, 
we are able to estimate fi a symp, the fractional density of matter that cannot be determined 
through gravity, to be 0.197-|-o.oi7; this is nearly equal to the measured value of the fractional 
density of nonbaryonic (dark) matter in the universe fl m — Qb 
WMAP [4j. Correspondingly, we estimate fioyn, the fractional density of matter in the 
universe that can be determined through gravity, to be O.tMl^Qg'j 1 ; this is nearly equal to 
the value of 0.0416±o;oo39 for Q B measured by WMAP [4j. We have also calculated a s , the 
rms fluctuation in the fractional density of matter within a distance of 8k" 1 Mpc, as a direct 
check of our model. Using the average rotational velocity and core sizes of 1393 galaxies 
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33j, 
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37| that 



obtained through four different set of observations 15, 
span 25 years, we obtain the value of 0.68±o.n for Cs; this value is within experimental error 
of the measured value of O^ll^g by WMAP. Finally, we have calculated -R200, the radius 
of the galaxy at which the density equals 200 times that of the critical density, to be 270±i 30 
kpc; this value also agrees well with observations. 

Interestingly, fi asymp /nA depends only on the dimensionality and symmetry of spacetime, 
and the exponent a\. This suggests that there is an underlying coupling between Dark 
Energy and matter in the theory. Such a coupling has been dismissed before, primarily 
because it is believed that the coupling would result in a "fifth force" that would already 
have been observed [8j . The results of this paper suggest that with a suitable choice of this 
coupling, its effects will not be currently measurable. 



A summary of the results here have appeared elsewhere [381 ] . Here, we provide the details 
of the theory and the calculations. 



II. EXTENDING THE GEODESIC LAGRANGIAN 



While there is no consensus as to the nature of Dark Energy — whether it is due to 



the cosmo logical constant A DE or to quintessence [39 
equations 

— c 



411 ] — modifications to Einstein's 
8%G 



-T 



(2) 



(where T^ v is the energy-momentum tensor for matter, R^ v is the Ricci tensor, R is the 
Ricci scalar, Greek indices run from to 3, and the signature of g^ v is (1,-1,— 1,-1)) to 
include the cosmological constant are well known, and minimal. We only require that Ade 
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changes so slowly that it can be considered a constant. We note also that the action for 
gravity+matter is a linear combination of the Hilbert action with a cosmological constant 
term, and the action for matter. Any change to the equations of motion for test particles 
can thus be accounted for in T^, and will not change the form of Eq. ([2]). 
With the geodesic Lagrangian 



dx^ dx 



1/2 



C = mc^ u ——) , (3) 
and the GEOM 

v v V v v» ee — = 0, (4) 

(where v v = x u is the four- velocity of the test particle), it is straightforward to see that 
in the absence of Dark Energy, Eq. (jlj) is the most general form that a second-order evo- 
lution equation for a test particle can take that still obeys the equivalence principle. Any 
extension of requires a dimensionless, scalar function of some a fundamental property of 
the spacetime folded in with some physical property of matter. In our homogeneous and 
isotropic universe, there are few opportunities to do this. A fundamental vector certainly 
does not exist in the spacetime, and while there is a scalar (the Ricci scalar R) and three 
tensors (# M „, the Riemann tensor R^ a/ 3, and the Ricci tensor R^), R^a/3 has units of 
inverse length squared. It is possible to construct a dimensionless scalar m 2 G 2 R/c i for the 
test particle, but augmenting £ using a function of this scalar would introduce additional 
forces that will depend on the mass of the test particle, and thus violate the uniqueness 
of free fall principle. It is also possible to construct the scalar g^v^v u /c 2 , but because of 
the mass-shell condition v^v^ 1 = c 2 , any such extension of Cq will not change the GEOM. 
Scalars may also be constructed from and powers of R p v taA p by contracting them with 
the appropriate number of t^/c's, but these scalars will once again have dimension of inverse 
length raised to some power, and, as with the Ricci scalar, once again the rest mass m is 
needed to construct a dimensionless quantity. 

For a nonzero Dark Energy, the situation changes dramatically. With a universal length 
scale Xde, it is now possible to construct from the Riemann tensor and its contractions 
dimensionless scalars of the form, 

c 2 R R^ u v^v v c 2 v IJ 'v u a0 v fl v u v 1 vs 



A DE G' A DE G ' (A DE G) 2 MQ,P7 " ' (A DE G) 
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While extensions to Cq can be constructed with any of these terms, we are primarily inter- 
ested in the nonrelativistic, linearized gravity limit. In this limit, the first two terms are 
equivalent to one another, while the other terms are smaller than the first two by powers 
of R, and can be neglected. We therefore focus solely on the first term, and arrive at the 
extension 

C Ext = mc yi + © ( Rc 2 /A DE G) \ (g^—^-j = X[Rc 2 /A DE G]C , (6) 

for Cq with the added constraint that v 2 = c 2 for massive test particles, and v 2 = for 
massless test particles. 

If D(x) is the constant function, then ^Ext differs from £ by an overall constant that can 
be absorbed through a reparametization of time. Only non-constant D(x) are relevant; it is 
how fast D(x) changes that will determine its effect on the equations of motion. Indeed, in 
extending Co we have essentially replaced the constant rest mass m of the test particle with 
a curvature-dependent rest mass m£K [Rc 2 /A EE G]. All dynamical effects of this extension 
can therefore be interpreted as the rest energy gained or lost by the test particle due to the 
local curvature of the spacetime. The scale of these effects is of the order of mc 2 /L, where L 
is some relevant length scale of the dynamics, and thus the additional forces from C^t are 
potentially very large. For these effects not to have already been seen is for D(Rc 2 / 'Ad E G) 
to change very slowly at current experimental limits. 

A. The Extended GEOM for Massive Test Particles 

For massive particles, the extended GEOM from C^y± is 

c 2 (g^ -^-)V u \ogX[Rc 2 /A DE G], (7) 



dt 2 V c 2 

,2 _ „2 



where we have explicitly used v = c . It has a canonical momentum with a 



2 a 2 2 

P = PpP = m c 



l+®(Rc 2 /A DE G) 



(8) 



and the interpretation of mdl[Rc 2 / Ad E G] as an effective rest mass can readily be seen. 

The dynamical implications of the new terms in Eq. (JTj), along with the conditions under 
which they are relevant, can most easily be seen after noting that R = 4A EE G / c 2 + 8nTG / c 4 , 
where T = T£. Then m[Rc 2 /A DE G] = SH[4 + 8nT/A DE c 2 }, where the '4' comes from the 



dimensionality of spacetime. It is readily clear that in regions of spacetime where either 
T^u = or when is a constant, the right hand side of Eq. ([7j) vanishes, and our extended 
GEOM reduces back to the GEOM. 

Beginning with Einstein 29J, there have been a number of attempts to show that the 



GEOM are a necessary consequence of the Einstein's equations Eq. (j2J). Modern attempts 



at demonstrating such a linkage [30, |3jj focuses on the energy-momentum tensor, and make 
the assumption that the strong energy condition holds: G^ v t^t' v < (for our signature for 
the metric), where G^ u is the Einstein tensor, and t M , t' u are two arbitrary, time-like vectors. 
They also assume that the background metric remains fixed during the passage of the test 
particle. With this assumption, the background metric decouples from the motion of the 
test particle, and can be treated separately. From the dependence of the extended GEOM 
on the energy-momentum tensor T — which includes a contribution from the motion of the 
test particle itself — it is clear that this assumption does not encompass our extension of 



the GEOM. It is thus not precluded by [30|, l3jj- Indeed, we will explicitly construct the 



energy-momentum tensor for dust within the framework of the extended GEOM in Section 
II.D. 

B. Dynamics of Massless Particles 

For a massless particle, the equations for motion from £Ext is 



v" 



V u (W[4 + MT / Kdec 2 ^) = (9) 



By reparametizing dt -> 9t[4 + 8vrT(x) / ' K DE c 2 }dt [4j|, Eq. © reduces to v u V u v fl = 0. 
With the correct choice of parametization, zero-mass particles still obey the GEOM. The 
usual general relativistic effects associated with photons — the gravitational redshift and the 



deflection of light — are thus not effected by our extension of the GEOM. This result is to 
be expected. Photons are conformal particles, and as such, do not have an inherent length 
scale to which effects can be compared [3]. 

C. Impact on the Equivalence Principles 



The statements 
following: 



431 ] of the equivalence principal we are concerned with here are the 
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Uniqueness of Free Fall: It is clear from Eq. ([Tj) that the worldline of a freely falling test 
particle under the extended GEOM does not depend on its composition or structure. 

The Weak Equivalence Principle: Our extension also satisfies the weak equivalence principle 
to the same level of approximation as the GEOM. The weak equivalence principle is based 
on the ability to choose a frame near the worldline of the test particle where r^ Aj g « 0; the 
Minkowski metric, r]^, is thus a good approximation to g^ v in the neighborhood around it. 
However, as one deviates from this world line corrections to r]^ u appear, and since a specific 
coordinate system has been chosen, they appear as powers of the Riemann tensor (or its 



contractions) and its derivatives (see 43| and 44j). This means that the larger the curvature, 



the smaller the neighborhood about the world line where r]^ is a good approximation of 
the metric. Consequently, the weak equivalence principle holds up to terms first order in 
the curvature, and since the additional terms in Eq. (J7J) are linear in R, our extension of the 
GEOM satisfies the weak equivalence principle to the same order of approximation as the 
GEOM does. 

The Strong Equivalence Principle: Because we only change the geodesic Lagrangian, all 
nongravitational forces in our theory will have the same form as their special relativistic 
counterparts. 

D. The Energy-Momentum Tensor 

As we have changed the equations of motion of test particles, we would expect the energy- 
momentum tensor for test particles to change as well. To see how it changes, we begin with 
the usual tensor for an inviscid fluid with density p, pressure p, and fluid velocity v^x): 

Tfj, u = pv^v v - (g^ - Vj ^~) V- (10) 
This form for T„„ depends only on the spatial isotropy of the fluid, and holds for both 



the GEOM and the extended GEOM. Following 43J, energy and momentum conservation, 
VT^j, = 0, requires that 

= v„V(p + p/c 2 )v^ + (p + p/c 2 )V u v v v^ + (p + p/c 2 )v v V v v^ - V„p. (11) 
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Since v 2 = constant even within the extended GEOM formulation, projecting the above 
along gives once again the first law of thermodynamics 

diVpc 2 ) = -pdV, (12) 

where V is the volume of the collection of particles. As such, the standard analysis of the 
evolution of the universe follows much in the same way as before under the extended GEOM. 

Next, projecting Eq. (ITT]) long the subspace perpendicular to t> M , gives the relativistic 
version of Euler's equation 

= (p + |) v»V„v, - (g» v - V> (13) 

We are concerned with the motion of matter in galaxies, and for such a system, test particles 
do not interact with one another except under gravity. This corresponds to the case of dust. 
If test particles in the dust follow the GEOM, then from Eq. ([TBI . T^°" Dust = pv ll v u and 
p = 0. On the other hand, if the test particle follow the extended GEOM, the situation 
changes. Using Eq. ([7]), Eq. ( fl~3l) becomes 

= (g, v - ^) { (pc 2 + p) V log m - V>} , (14) 

so that 

(pc 2 + p) V M <K - <KV M p = i^DECv^ (15) 

where £ is a constant. By contracting the above with v^, it is straightforward to see that 
if £ 0, p will increase linearly with the proper time. This would be unphysical, and we 
conclude that £ must be zero. 

As we are interested in the nonrelativistic, linearized gravity limit, T = pc 2 — 3p ~ pc 2 ; 

is a function of p only in this limit, and so, consequently, is p. Equation ( fTBT) then results 

in 

M = "V + c 2 *[4 + 8,P/A DE ] I m + Sns/AoE y (16) 

Given the density, the pressure — and thus the energy-momentum tensor for dust, T Ext ~ Dust , 
under the extended GEOM — is determined. 

To determine p, we note from Eq. (1161) that p ~ p 2 for p — > 0, while p ~ Ade when p ^> 
A.de/87t. We may thus still approximate T®*^^ « pVpV v in the nonrelativistic, linearized 
gravity limit. We next perturb off the Newtonian metric r]^ u through = + h^ u , where 
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the only nonzero component of h^ v is Hqq — 2$/c 2 , and $ is the Newtonian potential. It 
satisfies 



V 2 $ + 2- 



A EE G 



$ = AnpG - A DE G, 



(17) 



in the presence of a cosmological constant. 

As usual, the temporal coordinate, x°, for the extended GEOM is this limit is approxi- 
mated by ct to lowest order in |v|/c. The spatial coordinates, x, on the other hand, reduces 



to 



fx 
~a¥ 



-V$- 



Airc 2 



2) / (4 + 8vrp/A i?E ) 



Vp. 



(18) 



A DE J [1 +2) (A + 8tt P /A de )_ 
In principle, p can then be determined through the collection motion of the stars within 
galaxies. 



E. A Form for ®(x) and Experimental Bounds on a\ 

Since our extension of the GEOM does not change the equations of motion for massless 
test particles, we expect Eq. (J7J) to reduce to the GEOM in the ultrarelativistic limit. It 
is only in the nonrelativistic limit where deviations from geodesic due to the additional 
terms in Eq. ([7]) can be seen. We therefore focus on the impact of the extension in the 
nonrelativistic, linearized-gravity limit of Eq. ([7j), and begin by constructing D(x). 

For the addition terms from the extended GEOM not to contribute significantly to 
Newtonian gravity under current experimental conditions, S}'(4 + 8ttp/Ad E ) —* when 
p >> Ade/2tt. Note also that in the absence of the additional terms the motion of stars in 
galaxies is governed by a Newtonian, 1/r potential; what is instead observed is a weaker, 
logarithmic potential. These additional terms in the extended GEOM should thus contribute 
to the equations of motion of a test particle as though they were from a repulsive potential; 
this requires D'(x) < 0. 

The simplest form for D'(x) with these requirements is 

®>( x ) = — (19) 

where x is a normalization constant 

i r°° ds 



X J 1 + s x +«a 



(20) 
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To prevent negative effective masses, must be positive, so that 

^)=«f ^ (21) 

where a a > for the integral to be defined. While the precise form of D(x) is calculable, 
we will not need it. Instead, because 8np/ A DE > 0, 

V(4 + 8np/A DE )= X J2 7TZ L 4 + ^ ' < 22 ) 

^n(l + a A )+a A \ A DE ) 

?1=0 X 

while 

1 = 1 + 2^ ( -1 )" f23l 



Notice that in the a a — ► oo limit, D(x) — > 0, £s x t — ► £o, and the GEOM is recovered. 

Bounds on a a will be found below. For now, we note that for a\ > 1, X ~ 1 an d 
2) (4 + 87cp/A DE ) « 0. Thus, 

From WMAP, A^s = 7.21^q 2 84 x 10 -30 g/cm 3 , which for hydrogen atoms corresponds to a 
number density of ~ 4 atoms/m 3 . As such, the density of both solids and liquids far exceed 
A DE , and in such media Eq. (l2~4"j) reduces to what one expects for Newtonian gravity. Only 
very rare gases, in correspondingly hard vacuums, can have a density that is so small that 
the additional terms are relevant. To see when this may occur, consider the hardest vacuum 
that we know of at ~ 10 -13 torr 46) . For a gas of He4 atoms at 3 deg K, this corresponds to 



a density of pn m it ~ 1CT 18 g/cm 3 . Even though pn m it is still 11 orders of magnitude smaller 
than A DE , because the scale of the acceleration from the additional terms in Eq. ( [24]) is so 
large, effects at these densities can nevertheless be relevant. 

Let us consider an experiment that looks for signatures of the extension of the GEOM 
Eq. (!24j) by looking for anomalous accelerations (through pressure fluctuations) in a gas 
of He 4 atoms at 3 deg K, and p = p y m \ t- Inside this gas we consider a sound wave with 
amplitude epn m it propagating with a wavenumber k. Suppose that the smallest measurable 
acceleration for a test particle in this gas is abound- Then, for the additional terms in Eq. ( 1241) 
to be undetectable, 

^c 2 xf A DE V A 

abound > — „ ek. (25) 

2 V 87rpi imit / 
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This gives a lower bound on «a as 

_ log[2a bound /c 2 xeA;] 

^Abound — "j FT 73 f" • 1^0 

log [A DE /8irpi imit \ 

For e = 0.1, k = 1 cm -1 , and abound = 1 cm/s 2 , aA bound ranges from 1.28 for Ade = 10~ 32 
g/cm 3 to 1.58 for A DE = 10~ 29 g/cm 3 . 

In idealized situations such as Einstein's analysis of the advancement of perihelion of 
Mercury, the energy-momentum is taken to be zero outside of a massive body such as the 
Sun; the right hand term in Eq. (JTj) will not clearly not affect these analyses. This argument 
would seem to hold for all other experimental tests of general relativity as well. It is an 
argument that is too simplistic, however. In practice, the T^ v in each of these tests does 
not, in fact, vanish; there is always a background density present. Except for experiments 
involving electromagnetic waves, what is needed instead is a comparison of the background 
density with Aoe- It is only when this density is much greater than Aoe/^ that the 
additional terms in Eq. (JTJ) will be negligible. 

We have seen when this condition for the density holds for terrestrial experiments. Con- 
sidering now the traditional tests of general relativity, only in experiments involving motion 
of massive particles — such as the motion of Mer cury or the state-of-the art Eovtos-type 



experiments done recently by Adelberger 47|, |48|, |49] — will the effects of the extension be 
seen. The number density of matter at Mercury's orbit is roughly 100 atoms/cm 3 , however, 
corresponding to a mass density greater than ~ 10~ 23 g/cm 3 , which is orders of magnitude 
greater than A DE . It also has a corresponding length scale, c/ y/pG, that is on the order of 12 
Mpc, which is orders of magnitude larger than the size of the Solar System. The additional 
terms in the extended GEOM thus cannot appreciably affect the motion of Mercury, or 
any other solar body. Next, while the pressures under which Adelberger 's experiments were 
performed were not explicitly stated, as far as we know these experiments were not done at 
pressures lower than 10~ 13 torr; we would not expect effects from the additional terms to be 
apparent in these experiments either. We therefore would not expect the effects of Eq. (T53|) 
to have already been seen experimentally. Instead, with the average galactic-core density 
~ 10~ 24 -10~ 22 g/cm -3 and sizes of galaxies ~ 100 kpc, it is on the galactic length scales 
and longer where our extension to the GEOM will become important, and its effects felt. 
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F. Connections with Other Theories 



As unusual as the extended GEOM Eq. © may appear to be, there are connections 
between it and other theories. 

1. The Class of Scalar Field Theories in Curved Spacetimes 

The Klein-Gordon equation corresponding to the extended GEOM is 

m 2 c 2 ( Rc 2 \ 

W + ^ (l + 2(4) - 4JD'(4) + S'(4)£^J = 0. (27) 

where we have expanded D(i?c 2 / ' K DE G) about R = 4A DE G/c 2 . Although the relativis- 
ts Klein-Gordon equation for a scalar field theory can be straightforwardly generalized to 
curved spacetimes, it has also been generalized as 



m 2 c 2 



V 2 <p R +[^ r + ZR) ( j )R = 0, (28) 



since for £ = 1/6 the scalar field will be conformally invariant even though m/0 [451 ] . The 
similarity between Eqs. (1271) and (1281) is readily apparent. 

2. MOND 

As with MOND, the addition terms in our extension of the GEOM are nonzero at galactic 
length scales, while on terrestrial or interplanetary scales they are negligibly small. Like 
MOND, our extension is able to explain the galactic rotation curves, as we show in the next 
section. Within the MOND theory there is a fundamental acceleration scale cxmond = cH 



261 ]: in our analysis the scale that measures the additional contributions to the GEOM 



in Eq. (l24"j) is a Ext ~ c 2 /L. As galactic rotation curves are driven purely by gravitational 



effects, X EE is the only natural length scale, and thus aExt = c^GA DE ; numerically agxt ~ 
omond- Our extension of the GEOM thus gives an explanation for both the modification 
of Newtonian gravity that MOND proposes, and the fundamental acceleration scale that 
appears in theory. 

However, unlike MOND our extension of the GEOM is done within the framework of 
general relativity, and still requires the existence of Dark Matter. Although at the nonrela- 
tivistic, Newtonian level there is no separation between the force of gravity and the response 
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of matter to it, in general relativity there is. Our theory, in keeping the form of Eq. (j2J), does 
not change how matter curves spacetime; it changes how spacetime affects matter. Massless 
test particles still travel along null-geodesics, and obey the GEOM; the gravitation redshift, 
the deflection of light by massive objects, and gravitational lensing are all not affected by our 
extension. This is not the case for MOND, which was proposed at the Newtonian gravity 
level as a new theory of gravity. The theory must not only be extended to a relativistic 
one, but the response of electromagnetic fields must be extended as well, and this extension 
must be done in such a way that effects such as the red shift and gravitational lensing are 
unchanged. This program is not needed in our approach. 



3. The f(R) Theory 



Proposals for introducing additional terms of the form f(R) to gravity has been made 



before (see 



50j and 5J| for reviews), but at the level of the Hilbert action for R. These 



theories were first introduced to explain cosmic acceleration without the need for Dark 
Energy 



52 



531 ] using a 1/R action, and further extension have been made [54l . |55j . They are 



now being studied in their own right, and various functional forms for f(R) are now beinj 
considered. Indeed, connection to MOND has been made for logarithmic f{R) terms 561.157 [ 



59 



60 



61 



and with other choices of f(R), connection with quintessence has been made [58 
621 ] as well. Importantly, issues with the introduction of a "fifth force", and compatibility 
with terrestrial experiments have begun to be addressed through the Chameleon Effect (see 



64 



65. 



66, 



67| and an overview in }58|). This effect is a mechanism for hiding the effects of 



field with a small mass that would otherwise be seen. 



III. DARK ENERGY AND GALACTIC STRUCTURE 

While definitive, a first principles calculation of the galactic rotation curves using Eq. ffTB]) 
to describe the motion of each star in a galaxy would be analytically intractable. Instead, 
the approach we will take is to show that given a model of a stationary galaxy with a specific 
rotation curve, we are able to derive the mass density profile of the galaxy. The logarithmic 
interaction potential observed for the motion of stars in the galactic disk follows. We then 
will show that an idealized pseudoisothermal density profile will result in a lower free-energy 
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state for the calculated density than an idealized, cuspy-profile density. 



A. A Model Galaxy 



A number of geometries have been used to model the formation of galaxies (see 12|). 
Because we will be making connection with cosmology, we are interested in the large-distance 
properties of the density profile, however, and at such distances the detail structures of 
galaxies are washed out; only the spherically symmetric features of the galaxy survive. We 
thus use a spherical geometry to model our idealized galaxy, and divide space into the 
following three regions. Region 1 = {r \ r < r^, and p ^> Ae>e/2tc}, where r# is the galactic 
core radius. Region 11 — {r\r > r#, r < rjj, and p ^> Aoe/2tt} is the region outside the core 
containing stars undergoing rotations with constant velocity; it extends out to a distance of 
rn, which is determined by the theory. A Region III = {r | r > rn, and p <C Ade/^} 
naturally appears in the theory as well. 

We assume that all the stars in the model galaxy undergo circular motion. While this 
is an approximation, galactic rotation curves are determined with stars that undergo such 
motion, and we use these curves as inputs for our analysis. The acceleration of each star, 
a = x, is then a function of the location, x, of the star only. As such, we can take the 
divergence of Eq. (fT8|) . and obtain 

w-'-^'-j^gGeM- (29) 



where 

- 2 (p) = 7-^1+(4 + ^) h (30) 



and /(x) = —V • a/47rG is considered to be a driving term. Because we are dealing with 
only gravitational forces, we do not differentiate between baryonic matter and Dark Matter 
in p. 

In deriving Eq. (1291) . we used the Newtonian relation V 2 $ = 4irpG instead of the full 
expression Eq. (TTTj) . We do so because A DE is so small that it may be neglected for most 
of the regions we are interested in. While the contribution 2AdeGQ/c 2 to Eq. 017|) means 
that <l> oscillates, it does so on a length scale \de/\/2, which is longer than \ H /\/2; we will 
find that p is exponentially small where this scale is relevant. We are also mostly interested 
in regions where p » A DE /27i, and in this region the term —A DE on the right-hand-side of 
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Eq. (ITT)) is negligible; where p ~ Kde/2ti is precisely where p — > exponentially fast. The 
A^E term in Eq. (ITT)) therefore only insures that $ — > — c 2 /2 asr — > oo, and this can be 
taken into account with the appropriate boundary conditions for $. 
It is straightforward to see that 

V-a=-~[n;(r) 2 ], (31) 

where t>(r) is the velocity curve for the galaxy. Thus, given a v(r), a can be found and 
/(x) determined. For the observed velocity curves, we use a particularly simple idealization: 
v ldeal (r) = vht/th for r < r#, while -y ldeal (r) = vh for r > r#, where vh is the asymptotic 
value of the rotation curve. While y ldeal (r) is continuous, f(r) is not, and we find that 
f(r) = 3v 2 H / A7tGr 2 H = pn for r < r# and f(r) = v 2 H j^Gr 2 = pHrfj/3r 2 for r > r#. 

Analytically, this idealized velocity curve t> ldeal (r) is more tractable then the velocity curve 



for the pseudoisothermal density profile profile 



v^ir) = JAttGphR 2 



c 



1- Cretan, 

r \R C 



(32) 



However, because it has the same limiting forms in both the r <C r# and r ^> r# limits, 
v ldeal (r) functions an idealization of t)P _lso (r) as well. For the r 3> m limit, we need only 
identify y 1 \-nGpnR 2 c = and in the r -C r# limit we need only identify r# = a/3-Rc- 

For density profiles with a cusp-like structure Eq. OH), the situation is more complicated. 
Here, it is the density profile that is given, even though it is the velocity curve that is 
observed. While it is possible to integrate the density profile Eq. (TTJ to find the corresponding 
velocity curves f CU sp( r ), both the maximum value of v cusp (r) and the point where its slope 
changes — giving the size of the core — are different depending on the density profile used. 
Without a value of the core size that is consistent from one cusp profile to another, it is not 
possible to compare profiles and their free energies. 

Although it is possible in principal to determine the core sizes for each of the cusp profiles, 
doing so will be analytically intractable. Instead, we account for the different density profiles 
by taking 

/(x) = ( (r "' r)1 ifr - rH ' (33) 
[ \ph {th/tY if r > r H . 

The core size is set to be r H , and for the specific case 7 = and (3 = 2, Eq. (1331) reduces 

to the /(x) for v ldeal (r). For the velocity to be finite at r = 0, 7 < 2, while for it be finite 
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as r — > oo, (3 > 2. We will see that the density profile calculated from /(x) goes as 1/r 7 
for r < ru and as for r > r#; it is necessarily continuous at r = r#. Both limiting 
behaviors are equivalent to Eq. (pQ) , and thus Eq. (|3"3l results in an idealization of the density 
profiles in Eq. (pQ), while allowing us to consistently compare profiles. 



B. Observational Bounds on cka 

An estimate for a a can be obtained by comparing two length scales. From Eq. (129j) . the 
density near r# decreases with characteristic length scale 

A corc - , ■ , = X 1/2 \oe ^ . (34) 

K(p H /A DE ) \8irp H J 

The size of the galactic core r H should be proportional to A core , since at distances much 
smaller than th galactic dynamics are driven by Newtonian gravity, while at distances much 
larger than r# they are driven by the extended terms in Eq. (1241) . Fixing the A core in Eq. (1341) . 
we obtain an estimate for 

= log(xAWALJ _ 

log(87rp H /A DS ) 1 ; 

Figure [1] shows graphs of «a in term of A DE with the core density fixed at p H = 1CT 24 
g/cm 3 , and various values of A core . The characteristic length A core cannot exceed Th, nor can 
it be too much smaller than it; the values of A cor e chosen in Fig. 1 reflects this. Graphed also 
is the lower bound on set in section II. D. This bound, combined with Eq. (1351) . brackets 
«a within the triangle bound by 32 pc < A core < 100 pc and A DE > 3.3 x 10~ 31 g/cm 3 , and 
limits ctA > 1.35; A E e — 7.21 x 10~ 30 g/cm 3 lies within this triangle. Given this result, we 
take «a — 3/2 as the representative value of a>\ for this part of the paper. A definitive value 
for a\ will be set in Section V. 



IV. THE DENSITY PROFILE OF THE MODEL GALAXY 

In both Regions I and II, pu ^> Aoe/^ and Eq. (1291) may be approximated as the 
following, 

/(») = P - (^) {v^P - (1 + "a) 1 ^} . (36) 



18 




FIG. 1: Graphs of a a with respect to Ade at pu = 10 24 g/cm 3 for various values of A CO re < ru- 
The graph of the lower experimental bound Eq. (|26p is also included. 

where u = r /x 1 ^ 2 ^de, and V u denotes derivative write respect to u. The solution to Eq. (13S|) 
minimizes the following functional of the density: 

2 



1 


vj 




2a A 






A DE 


\ OA 


8nf(u) 


8np 




A-DE 



OA 



"A 



"A - 1 



A 



DE 



inp 



<*a-i 



+ 



(37) 



which we view as an effective free energy for the system. Here, we have chosen the scale as 
A-dec 2 . 

In Region III, on the other hand, the following linearization of Eq. (|29|) is appropriate, 

1 



= p- 



(38) 



1 -\- 41+o-a 

We will find from a detailed analysis of the solution for p in Region II that the driving term 
/(r) is negligibly small in Region III. A calculation of the free energy in this region will not 
be necessary. 
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A. The Solution for Region I 



For the <y ldeal (r) curve — corresponding to 7 = — it is clear that the only solution, pi(r), 
for Eq. (J3HD i n Region I is the constant solution pi(r) = p#. The free energy for this solution 
is easily calculated 



1 , A T" (^) < 39 > 

a A -l 6 \8irp H J 
For a general 7 > 0, perturbation theory is used to find solutions Eq. ( 1361) . We first scale 
rj = p/ph and j: = u/uh < 1, so that 



QA-l 



P 7 = tj + — Vtj 



where the small parameter 



1 fA 



1+«A 



~ 10- 2 -10" 9 , 



(40) 



(41) 



for a = 3/2, p# ~ 10~ 24 -10 -22 g/cm 3 , and r# ~ 1 kpc -8 kpc. There are two approaches 
to solving Eq. (|40p perturbatively. The first treats the V 2 term as a perturbation on the 
solution pi = f . Doing so gives 



(ft) / \ IT 

Pi r = Ph , 



1 - £7(1 + 7a A ) 



7(l+a A )-2' 



(42) 



It is valid when 7 > 2/(1 + a A ). For the second, we take 3 = e/rj^, so that 



? -7 = e i/«A 3 -i/a A + _V 2 j. 

"A 



(43) 



Treating e 1 / QA 3 1 / aA now as the perturbation, the solution for 3 that is finite at r. = gives 
- £ (2-7)(3- 7 ) xl/aA 



pf\r) 



-7 



i + i 



r 

V«a ( 2 - 7 )(3- 7 ) / r \ 2 ~^ 1+a ^ 



e(2- 7 )(3- 7 ) 



(44) 



(2-7)(3-7) W. 

where 7 = (2 — r y)/a\. The solution is now valid for 7 < 2/(1 +a\), which again corresponds 
to 7 >2/(l + « A ). 

Calculating the free energy for these two perturbative solutions follows straightforwardly. 
To lowest order in e, the free energy for pj is 



I T (a) 
- r 7>0 



1 A DE c 2 r 3 H ( A 



«a — 1 6 



8ttph 



Q A -1 



l + 7 (a A -l)/3' 



(45) 
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while the free energy for pj is 

i T b aAc 2 r 3 H { 8np H u H \ 2 

^>°- (2 - 7 ) (3 - 7 )4 *de ) ■ W 

It is clear that I J 7 1 =o < 1 ^>2/(i+a A ) < < 7 -^<2/(i+a A )' ^ ne idealized, pseudoisothermal 
profile corresponding to 7 = is thus the state of lowest free energy in Region I. Physically, 
this results because of the curvature term ~ |Vp| 2 > in Eq. ( 1371) . Just like the lowest free 
energy state of a Landau-Ginzberg free energy functional, this term only vanishes for the 
constant solution; for all other solutions it contributes positively to the free energy. 



B. The Solution for Region II 

In this region, 



-Ph — =Pii-X*de o Vp w -(l + a A ) ^, (47) 



1+Q A (- ivy^ |2 

„.., , s , {l+ajFez. 

3 v r J \8np n J { p n 

and we denote pn as the solution for p in Region II. We undertake an asymptotic analysis 
68| of Eq. ( 14"TI) by making the anzatz that within Region II there exists a point r asymp 

asymp ; 



beyond which p_ff(r^/r)^/3 -C p(r). For r > r asymp , we can then neglect the driving term in 



Eq. (147|) . leaving the homogeneous equation 

= Pasymp - f jr^-) { V^ asymp - (1 + c A ) |VttPasymp|2 | . (48) 

1. Asymptotic Analysis and the Background Density 

We look for a power-law solution to Eq. fj48l with the form 

Pasymp = ^^S(a A )u P , (49) 

and find that 

p(l-q A p) 1 

U - 1 + - X l+a AM p(l+ QA )+2- ,0U) 

£(a A ) 



This gives p = —2/ (1 + ct A ), with £(a A ) the solution of 

2(l + 3a A ) 
;i + a A ) 2 [S(« A )] 
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Positivity of p a s ym p requires that £(«a) > 0, and this requires that there are positive solutions 
to Eq. flSTj) . Such solutions exist only if a A is the ratio of a odd integer with an even integer. 
Our choice of «a = 3/2 satisfies this criteria, and we arrive at the following asymptotic 
solution 



A 



P 



DE 



asymp 



8tt 



XA 2 



DE 



l/(l+a A ) 



where 



E(a A ) 



' 2(l + 3a A ) l 1/(1+QA) 
,(1 + «a) 2 . 

To justify our anzatz that r asymp lies within Region II, we set /(r asymp ) 
and find 



(52) 



(53) 



asymp y asymp J > 



asymp 
TH 



< 



o 2/(l+a A )\ ( 1+Q A)/2a A 



3A 



DE 



E(a A ) 



(54) 



where u H = r H /x 1,2 ^DE- For a A = 3/2, p H ~ 10 24 g/cm 3 , and r H = 1 kpc, r asymp < 
1.77r#; the anzatz is valid through the great majority of Region II for the range of galaxies 
we are interested in. The upper limit ru to Region II, on the other hand, is found by setting 
87rp asymp (r / /)/A jDi; = 4, which gives 

(l+<* A )/2 



ru 



X 



V2 A 



DE- 



(55) 



For a A = 3/2, r n « 0.20 A DS . 



T/ie ./Year Core Density 



Structural details of the galaxy cannot be seen from p asymp . Instead, we take pu 
P aS ym P + p\i an d expand Eq. fT4T|) to first order in 



2(l + 3a A ) 
(1 + «a) 2 



Ph_ (u h \0 
3 V u J 



2(1 + 3a A ) p/ 7 



+ v 2 pA. 



(56) 



(l + a A ) 2 « 2 

where pjj = u 2 p}j. In the special case [3 = 2, the particular solution to Eq. (15"6"1) is again the 
constant solution, but now for this corresponds to p}/_ / 3 =2 = P-ff r n /3r 2 , as expected for 
an idealized pseudoisothermal profile. 

Boundary conditions for p n are se t & t the r = th surface, 



Ptf = P. 



asymp 



{ru) 



PlIVH) 







dp, 



asymp 



H 



du 



3_ 

du V u 2 



Ph 



(57) 
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where we have made use of the result of the Region I free energy analysis and set 7 
The solution to Eq. fl56l) for these boundary conditions is 



COS 



VQ log 



+ C s i n sin 



i/ log — 



r 



0. 



(58) 



where v = [2(1 + 3a A )/(l + « A ) 2 - 1/4] 1/2 , and 

»l + 1/4 



.4, 



ul + (5/2 - py 



Ph 
5 



I a - A DjE S(q a ) 
3 ppH 8tt ^2/(i+«a) 



r ^p H (5/2-/3) 



1(1 + 5a A ) A DE S(q!a) 



(59) 



3"^ riJV " / 2 (1 + « A ) 8tt Mh 2/(i+« a )- 

The density pn(r) thus consists of the sum of two parts. The first part, p a s y mp( r ), 
corresponds to the background density, and depends solely on Dark Energy, fundamental 
constants, the exponent qa, and the dimensionality and underlying spatial symmetry of the 
spacetime. It is universal, and has the same form irrespective of the detailed structure of the 
galaxy. The second part, p}j(r), does depend on the detail structure of the galaxy. Variation 
in p]j are measured on a scale set by Th, the core size, in contrast to p a sym P , whose variations 
are measured on a scale set by A^, the Dark Energy length scale. While our analysis is 
done only to first order in the perturbation of pu, this feature of p x n holds to higher orders 
as well. 

The perturbation, p}j, itself depends on two terms. The first has a power law dependence 
of (th I r)^, while the second has a power law dependence of (r# / r) 5 / 2 . Thus, near the galactic 
core p]j ~ r~ q , where q = max(/3, 5/2); to this level of approximation, the density profile 
near the core varies at least as fast as 1/r 5 / 2 , irrespective of whether the density profile is 
cuspy or pseudoisothermal. For the r 3> r H the opposite is true, and now q = min(/3, 5/2); 
pjj decreases no faster than 1/r 5 / 2 . 
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3. Free Energy Analysis 



The free energy for the density pn separates into three terms: 11 T = II T asyrav + 
// ^ r asym P -/3 + 11 , where 

A DE c 2 



asymp 



8. ^L/^ 



A 



DE 



8np 



asymp 



OA 



asymp 



4alA DE c 2 (x 1/2 \DE) 3/2 F(^a)] 2 M (i+5a A )/(i+q A ) (6Q) 
" K-l)(l + 5a A )[S(a A )] 2(1+aA)M// 
is the contribution to the free energy due to the background density only. The integration 
is over Region II, and uu = i"ii/x 1 ^ 2 ^de- 
The second term is 



n 



T 

j f 



asymp— /3 



( X V2 A 



DE) 



I d*uf(u)(-^L_ 

<D n \ 87r Pasymp 

„2 



+ 



87TQ;aC 



A DE [S(a A )] 2 d+°A) . 



a A (8v 



2 2 
i C 



A DS (l + a A )[S(a A )] 2 d+«A) 
A DS c 2 (l + a A ) 2 / Anp H 



u p // (M)Vp asymp • dS. 



{u 5 H p] I {u H )p 

asymp 

(u H ) -u] I p) I (u II )p 

asymp 



£(a A )i4 



5-2/(l+a A )-/3 
7/ 



5-2/(l+a A -/3) ' 



,(61) 



2(l + 3a A ) V3A DS ; v y a ^ 5-2/(l + a A ) ~P 
where dD n is the boundary of D n aX r = r H and r = r n . This is the contribution to the 
free energy due to the interaction between p a sym P and /(r). It is straightforward to see that 

-{uh/ujjY if /? < 5/2, 

-K/M 5/2 if 5/2 < (3 < 5 - 2/(1 + a A ). (62) 
±u 5-2/(l+a A ) if 5 _ 2/(1 + « A ) < ^ 

where the sign of the last term depends on the values of f3, a A , uh and pn- The magnitude of 
this term is very small, however, and it is clear that // ^ r a symp-/3=o < // -7~asymp-/3<5-2/(i+a A ) < 



n 



T 

J r 



asymp 



-P ~ < 



II jr 



asymp-/3>5-2/(l+a A ) • 

The third term 

a A (l + a A )c 2 { X 1/2 XdeY 



ii T i 



T l = 



2[S(a A )] 2 ( 1 +«A) 



4 1 / \ ^Pasymp , 
U Pij(U) — - • GtS + 

dD n Pasymp 
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a{x 1/2 \de) 3 f8nc 



[S(a A )]2(i+"A) \A DE 



(l + 3« A )(p/ 7 ) 2 2(1 + 3^) , 
+ -7^ —JWPn 



+ a A ) 2 u 2 



(63) 
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is the contribution to the free energy due to p x u only. As it ~ (p)/) 2 , this term is very small 
compared to the other two terms that make up the free energy, and can be neglected. We 
only note that like it is the constant solution to the differential equation that gives the 
lowest value of II J rl , but now the equation is for 'p} 1 = u 2 p 1 II , not p 1 II . This once again 
corresponds to (5 — 2, and as such, to a p x n ~ 1/r 2 . 

The total free energy, 11 T, in this region is thus smaller for (3 = 2 than for (3 > 2. 
Combined with the calculation for I J r , we conclude that the pseuodoisothermal rotational 
velocity curve will result in a density profile that gives the lowest free energy, and is the 
preferred state of the system. Other rotational velocity curves will result in density profiles 
that have a higher free energy. We therefore take 7 = and (3 = 2 for the rest of this paper. 



C. The Solution for Region III 

The solution to Eq. ( 1551) follows using standard methods, with the boundary condition 



pn(r n ) = p ni {r n ). The length scale y/x/0- + ^ 1+cxa )^de ~ 0.15 X DE for a A = 3/2, while 
t u ~ 0.20X DE ] to a good approximation r u ~ a/x/(1 + 4 1+a *)\ DE . As y/x/0- + 4 1+QA )A DS 
is a scale set by the theory, we shall use this last expression for rjj from now on. Note also 
that Regions II and III overlap, and our approach of solving the nonlinear partial differential 
equation is self-consistent. 

The only solution that is spherically symmetric and finite at r — > oo is 



Pm{r) = ^£ E(aA) ^EH (1 + 4l+ « 4) J<i-*va + -> exp L _ j_ Jl+^E] , (64) 



8tt 

and in this region the density decreases exponentially fast. 




D. The Effective Potential 



Note that Eq. ( fl51) may be written in terms of an effective potential, 5J(x), as x = — V23 
where 

9J(x) = $(x) + c 2 log9 : l[4 + 87rp/A Di; ]. (65) 

It is thus not the gravitational potential $(x) that determines the dynamics, but QJ(x). To 
see the implications of this, we begin by calculating $(x), but only for Regions I and II. In 
Region III, r > rjj, and motion in this region is unphysical. 
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Integrating V $ = AnGp gives in Region I 

*< r > = 7* (66) 

where the c 2 term comes from requiring $(rjj) = — c 2 /2 instead of $(oo) = 0. That this is 
the usual expression for the Newtonian potential can be seen from pn = ?>v 2 H j ^Gr 2 H . For 
Region II, 

W 4« A (l + 3a A ) F/n; (l + 3« A ) 6 V^y 



c 2 2 

h 

2 H 



1 + 4 



'1 + "A 



|2 



;i + 3a A )J 4~ ^ v iW (l + 3a A ) H 
f , (l + ^ A ) 2 _ c 2 (1 + a A )(l + 2a A ) 2aA/(1+aA) | 

\ b ^(l + 3a A ) " 2" XL(aA) (l + 3a A ) J r ' (67) 

The 1/r terms in $ are expected from Newtonian gravity, and is due to the boundary 
conditions for $ at r = r#; this is true for the constant terms as well. The logarithmic term 
is due specifically to the driving term /(r), as expected. It is a long-range potential that 
extends out to rjj, and could potentially explain the non-Newtonian interaction between 
galaxies and galactic clusters in addition to the galactic rotation curves. The p]i{r)r 2 term 
is due to the perturbation in the background density, and contains terms ~ 1/r 1 / 2 . It is 
due to both the boundary terms in p\ I and the boundary conditions for $(r). Finally, the 
u 2a A /{i+a A ) t erm j n gq ( pTj) i s d ue to the background density p a sym P , with origins rooted in 
the rest mass of the test particle. Note also that this term is less than c 2 , verifying that the 
nonrelativistic limit still holds. 

The c 2 term in $(r) increases as r 6 / 5 for a a = 3/2, and would dominate the motion of test 
particles in the galaxy if the extended GEOM depended on $(x) instead of 5J(x). Instead, 
it and the r 2 p]j term are canceled by the additional density-dependent terms in Eq. fl65l) . 
To see this, in Regions I and II p >> Ade/27t, and expanding Eq. (1651) gives 



2 n \r H J 2 2« A \8irp H 
in Region I, while in Region II, 

where we used p}j/p a symp << 1 and Eq. (I5T1) . The last two terms of Eq. (1691) cancels the first 
two terms of Eq. fIBTj) . and 5J(r) is simply the remaining five terms in Eq. fIBTl) ; the effective 
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potential thus increases only logarithmically as r increases. This is what is expected for a 
potential that determines the motion of the stars in galactic rotation curves. 

That the motion of stars in the galaxy is determined by 9J(x) and not $(x) has far 
reaching implications. The c 2 term in $(r) comes from the background density p asy mp, and 
thus the majority of the mass in a stationary galaxy does not contribute to the motion of test 
particles in the galaxy. Rather, it is the near-core density p}j that contributes to 23(x), which 
results in the very long-range, logarithmic potential that is observed. As such, observations 
of the rotational velocity curve of a galaxy will be able to determine the perturbation on 
the background density, pjj, but not p asym p itself. Consequently, since p asy mp(^) 3> p l n (r) 
when r ^> r#, the majority of the mass in the universe cannot be seen with these methods. 
In particular, the motion of stars in galaxies can only be used to estimate p(r) — p asym p( r ); 
the matter in p asymp (r) is present, but cannot be "seen" in this way. 

This behavior is expected for a background density. In the traditional theory of structure 
formation, the perturbation is off the average density of matter of the universe. This density 
is usually taken to be a constant, and thus cannot affect the motion of stars within the 
galaxy. It also fits well with our interpretation that the extension of the GEOM is due to the 
replacement of the constant rest mass m with a curvature-dependent rest mass. While the 
rest mass contributes to the Newtonian gravitational potential energy for geodesic motion, 
it is a constant, and does not contribute to the dynamics of the particle. We see a similar 
effect here. Our p asymp is not a constant, however. It increases as r — > r H , and this is a 
feature expected of cold dark matter. Indeed, we will find below that i? 2 oo is determined 
primarily by p asym p- As such, in the absence of all other forces the majority of the mass 
outside of the galaxy cannot be observed through its dynamics. Other means would have to 
be used. 

V. COSMOLOGY 

In the previous sections, we have focused on analysing the structure of a single galaxy. 
In this section, we will extrapolation these results to the cosmological scale, and perform 
a cosmological check of our theory. That this extrapolation can be done is based on the 
following two observations 

First, recent measurements from WMAP and the Supernova Legacy Survey put VL K = 
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-0.011 ±0 .oi2, while WMAP and the HST key project set Vt K = -0.014 ±0 .oi7 Q- In both 
cases, measurements have shown that the universe is essentially flat, and WMAP's deter- 
mination ofh = 0.732to!o32 

was made under this assumption, cLS IS cL ge of the universe to a t 

13.731 

a 15 Gyr. As such, the largest distance between galaxies is cto = R(Q)\h, where 
depends on the details of how the universe evolves, and thus on its thermal history; from 
measurements of both h and to, ^(fi) = 1-O3±o.o5- 

Second, the density of matter of our model galaxy dies off exponentially fast at r n . The 
extent of the mass of matter in the galaxy is thus fundamentally limited to 2rjj. This size 
does not depend on the detailed structure of the galaxy near its core; it is inherent to the 
theory. Moreover, as we can express 

.3n A (l 8 +V«A)] Xh ^ 

where Q\ = 0.716±o.o55 141 * s the fractional density of Dark Energy in the universe, we find 
that for ctA = 3/2, m = 0.52 Xjj. Thus, although the value of a\ was set to 3/2 by the 
structure of the galaxy on a galactic scale, the density distribution for the galaxy naturally 
cuts off at a radius equal to half to Hubble scale, which is precisely what is expected from 
cosmology. 

To accomplish the extrapolation, we use the properties of a representative galaxy for the 
observed universe to construct our model galaxy. This representative galaxy could, in prin- 
cipal, be found by sectioning the observed universe into three-dimensional, non-overlapping 
cells centered on a galaxy; given the spatial inhomogeneity of the distribution of galaxies, 
these cells will not all be the same size. Through a survey of these galactic cells, a represen- 
tative galaxy, with some average rotational velocity v* H and core radius r* H , can be found, 
and used to set the parameters for our model galaxy. While such a survey has not yet been 
done, there exists in the literature a large repository of measurements of galactic rotation 



curves and core radii 115 



36, 



37| . Taken as a whole, these 1393 galaxies are reasonably 



random, and are likely representative of the observed universe at large. 

While we were able to estimate of a a = 3/2 by looking at the galactic structure, the 
accuracy of this estimate is unknown; comparison with experiment will thus not possible. 
We instead require that r u = &(fl)\ H /2, which gives a>\ as the solution of &(Q) 2 (1+4 1+Q!A ) = 
327rx(a;A)/3fiA; this sets a\ = 1.51±o.ii- 
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A. ag, i?200? an d Galactic Rotation Curves 



Linear fluctuations in the density are denned as S(r) = Sp/(p) (69J, where (p) is the 
spatially-averaged density of matter in a set volume. The rms fluctuation of 5 is measured 
through cr| = ((£(x) 2 ) 8 , where the subscript denotes a spatial average over D 8 , a sphere 
with radius 8h~ x Mpc. In calculating a 8 , we choose as our origin the center of the unit 
cell containing the representative galaxy mentioned above. One estimate of the size of a 
typical unit cell as 5/i _1 Mpc, the characteristic length scale for the galaxy-galaxy correlation 
function [ooj ] : as this is only slightly smaller the 8/i _1 , it is reasonable to consider the density 
from only a single galaxy within D s . 

In our theory, p{r) varies significantly across D 8 . We thus begin by calculating 



(p)s = (pH0(r H - r))s + (pasymp (r)6>(r - r H )) 8 + (p l H {r)6(r - r H ))i 

3A D £ / 1+«A \ x -2/(l + Q A ) n / S 

= — (jr^r KK D{us) 

where 9(x) is the step function, us = Sh^Mpc/x^^^DE, and 



The function 

l + 3a A 

D(u 8 ) = l-y 8 1+aA + 



3y|(l + a A ) : 
2(l + 3a A ) 



c 



2w -2a A /(l+a A ) 



E(a A )x(aA) V c 



1 + 3a\ 
1 + «a 



P + 



3y 8 (l + a\f 



2(1 + 3« A ) 
sin [i/ logy 8 ] + 



^oCsin — ~ C c , 

2u 



UoCsin ~ 2° cos 



(71) 



(72) 



cos [f Q log 




where y$ = uh/u$, and 
C, 



= -f- I 7/ 2«a/(1+"a) 



3 6 V 1 + «a 



(74) 



Information on the structure of the galaxy is contained in (. As £ ~ 5 x 10" 3 for v* H = 200 
km/s, D{u 8 ) ~ 1, and it is the background density p aS ymp that dominates (p)&, and not the 
detail structure of the galaxy. 

This is not the case for (5(x) 2 ) 8 , which involves the integration of (pn) 2 over D s . Because 
(p}i) 2 ~ V r4 > ^ is now the behavior of the density near the core that is relevant. Indeed, 
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we find 



-1 + 



1 + 3«a 



3D(u 8 ) 2 \ l + a A 



1 + a a 



3a A 



+ 2C 



1 + «A 
oca - 1 



(Q A -l)/(a A +l) 



+ 



3 n ^ \ Lr («a - 3) ~ 

2 [_ 2(1 + o>a) 

1/2 I p, (a a - 3) 



(a A -l)/(a A +l) 



[l + a A )y 8 



AC 

y& 

9 
4 

4 



^oCcO! 
"(A 

3 



l + 3a A )/(l + a A ) 2 
1 



2(1 + a A , 

(«a - 3) 
" 2(1 + a A ) 



a. 



sin[i/ log?/8] + 



cos[i/ logy 8 ] + 



:C r 



+ 



Csin + C 2 + 



1 



1 + ^ 

1 



Csln — Ccos + 2z/oC cos C S ; 



(75) 



3 (a A + 1)(3o;a 
10 

l + (l + 3a A )/(l + a A ) 5 



1 + 3a A 
1 + a A 



Sa A C 



(a A - l)(3a A 



1/8 



5 + 



81 



4(1 + vl, 



(76) 



where for Eq. (1761) we have kept the lowest order terms in ( and ?/§. The first term in 
Eq. ( 1761) is due to the background density p a sym P - It depends only on a a, and contributes a 
set amount of 0.141 to a| irrespective of the structure of the galaxy. The last term is due 
primarily to the 1/r 2 term in p} 7 , and is due to the rotation curves. This term contributes 
the largest amount to erf, and depends explicitly on details of the structure of the galaxy 
through v* H and r* H . 

While there have been a many studies of galactic rotation curves in the literature, our 
need is for both the rotational velocity and the core radius of galaxies. This requires both a 
measurement of of the velocity as a function of the distance from the center of the galaxy, 
and a fit of the data to some model of the velocity curve. To our knowledge, this analysis 
has been done in four places in the literature [761 ] . While each of the data sets were obtained 



with similar physical techniques, there are distinct differences in their selection of ga 



in the exact experimental techniques used, and in their fits to rotation curves (see 71[ for a 



axies, 



new method of deriving the rotation curves from HI data). In fact, the Hubble constant used 
by each are often different from one another, and from the value of 73.2 km/s/Mpc given by 
WMAP. The reader is referred to the specific papers for details on how these observations 
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were made. Here, we only note the following: 



de Blok et. al. Data Set: de Blok and coworkers made detailed measurements of 60 LSB 
galaxies I?]], and fits to v p ~ lso (r) were done for 30 of them [ijj]. Later, another set of 
measurements of 26 LSB galaxies were made by de Blok and Bosma 16j, of which 24 



are different from the 30 listed in 



15]. Both the data for the 30 original galaxies, and 



the 24 subsequent galaxies are used here. Although the authors used various models for 
determining the mass-to-light ratio in their measurements, we will use the data the comes 
from the minimum disk model, as this was the one model used for all of the galaxies in the 
set. 

While fits to the NFW density profile were made in Q], we are primarily concerned 
with the fits by the authors to the pseudoisothermal profile velocity curve Eq. f)32p . As de 
Blok and coworkers were chiefly concern was with the density parameter pu for the profile 
and Rc, the fits were made with these two parameters. Standard errors for both pn and 
Rc were calculated and given. Our concern is with the asymptotic value of vh, however, 



and as this value is given by yA^fpjir^, we have calculated vh and its standard error from 
the published values of pu and Rc in [15-] to determine v* H for this set. The authors used a 
value of 75 km/s/Mpc for the Hubble constant. 



CF Data Set: In 



36j, Courteau presented observations of the velocity curves for over 300 



northern Sb-Sc UGC galaxies, and determined r# for each by fitting the curves to three 
different velocity curves, one of which is similar to the velocity curve for the pseudoisothermal 
profile used by de Blok and coworkers, 

v vcA (r) = -v c arctan ( - ] . (77) 
vr \r t J 

Like f p_lso (r), t> vcA (r) can be approximated by the idealized velocity curve y ldeal (>) used in 
our analysis. In the limit r ^> r t , v vcA « vc, which sets v c = Vh- In the limit r < rt, 
^vcA ^ vc(2r/iTr t ), which sets r t = 2rn/-n. Although there are differences between t> vcA (r) 
and yP _1S0 (r) within the two limits, our analysis here is based on the idealized velocity curve 
and all that is needed is the relationship between r#, and r t or Rc- 

The second fit was to a velocity curve where the steepness of the transition from the hub 
and the asymptotic velocity curves could be taken into account as well. Because of the work 
by de Blok and coworkers in 15], our focus here will be on the velocity curve Eq. ([77)1 . and 
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not this curve. The second velocity curve has a limiting form in the r — > limit that only 
agrees with our idealized profile in one special case, and this case does not hold for all the 
galaxies analyzed by Courteau using this profile. 

The third fit was to Persic, Salucci, and Stel's Universal Rotation Curve (URC) [72] . 
While the URC asymptotically approaches a constant velocity, at small r the URC has a 
r o.66 behave ; which is different from t> ldeal ; v P~ mo ^ anc [ -y vcA ; all of which varies linearly in 
the small r limit. We therefore did not focus on this velocity curve here. 

Values for vc and r t for 351 galaxies was obtained through the VizieR service 



(http://vizier.u-strasbg.fr/viz-bin/VizieR). The great majority of the rotation curves were 
based on single observations of the galaxy; only 75 of these galaxies were measured multiple 
times, with the majority of these galaxies observed twice. The data set reposited at VizieR 
contained these multiple measurements. We have averaged the value of Vc and r t for the 
galaxy where there are multiple measurements of the same galaxy. The standard error in 
the repeated measurements of a single galaxy can be extremely large; this was recognized 
in 36j]. A value of 70 km/s/Mpc was used for the Hubble constant by the author. 
Mathewson et. al. Data Set: In 37j], a survey of the velocity curves of 1355 spiral galaxies 
in the southern sky was performed. Later, the rotation curves for these observations were 
derived in [3] after folding, deprojecting, and smoothing the Mathewson data. Each of these 
velocity curves are due to a single observation. Courteau performed a fit of Mathewson's 
observations to the velocity curve Eq. (1771) for 958 of the galaxies in |36|] using a Hubble 
constant of 70 km/s/Mpc. The results of this analysis is reposited in VizieR as well. 



Rubin et. al. Data Set: In the early 1980s, Rubin and coworkers 32|, |33|, |34], |35J presented 
a detailed study of the rotation curves of 16 Sa, 23 Sb, and 21 Sc galaxies. This was not a 
random sampling of Sa, Sb, and Sc galaxies. Rather, these galaxies were deliberately chosen 
to span a specified range of Sa, Sb, and Sc galaxies, and as stated in [34|, averaging values 
of properties of galaxies in this data set would have little meaning. These measurements 
can contribute to the combined data set of all four measurements, however, which is why 
we have included them in our analysis. While values for the core radius were not given, 
measurements of the rotational velocity as a function of the distance to the center of the 
galaxy were; we are able to fit the data to the same pseudoisothermal rotation curve Eq. (|32|) 
used by de Blok, et. al. Results of this fit is given in Appendix A.l. A Hubble constant 
of 50 km/s/Mpc was used by the authors. 
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Wanting to be as unbiased and as inclusive as possible, we have deliberately not culled 
through the data sets to select the cleanest of the rotation curves. Nevertheless, we have had 
to removed the data for 27 galaxies from our analysis. A list of these galaxies and the reason 
why they were removed are given in Appendix B, where we have listed other peculiarities 
found with the data sets as well. 

The Vjj and r* H , and standard error for each, were calculated for each of the four data 
sets considered here. While Vjj is easily identified for all four, determination of r# is more 
complicated. For the de Blok et. al. data set, published values of Rc was first scaled by 
75/73.2 to account for differences in the Hubble constant used by the authors, and the current 
value of 73.2 km/s/Mpc measured by [4|. Then r H is obtained by using r H = V3R C . A 
similar calculation was made using the calculated values of Rc from Appendix A.l, but 
using 50/73.2 instead of 75/73.2 to account for differences in Hubble constants. For the CF 
and Mathewson et. al. data sets, published values of r t are first scaled by 70/73.2 to account 
for differences in Hubble constants, and r# is now obtained using r# = ixr t /2. A fifth data 
set is then constructed by combining the data from these four data sets. For each data set, 
Vh and r* H are then used in Eq. ( I75l) to calculate <7g; numerical derivatives of <7g were then 
used to calculate its standard error. Not surprisingly, A<7 8 is dominated by the standard 
error in a\. 

Results of these calculations are giving in Tabled! along with the t-test comparison of the 
calculated <Tg and Aa 8 with the value 0.761^Q;Q4g from [4J. Surprising, four of the five data 
sets give a value for erg that is within two-sigma of the WMAP value; they thus agree with 
the WMAP value at the 95% confidence level. The only data set that differs significantly 
from the WMAP value is the Rubin et. al. set, and it is known that for this set the selected 
galaxies are not representative of Sa, Sb. and Sc galaxies; this disagreement is thus not 
surprising. 

While the URC has a different power-law behavior at small r than t> ldeal ; v p~ iso ^ or -y vcA , 
the difference is small enough that it is unknown how <jg will change if the URC is used in 
its calculation instead of the t> ldeal used here. We leave this for future research. 

Given v* H and r* H , it is possible to calculate i? 2 oo by setting w 8 — > u 2 oo = R2oo/x 1 ^ 2 ^de in 
Eq. (ITT]) , and using v* H = 172.1 km/s and r* H = 11.82 kpc from the Combine data set. We 
numerically solved for this radius and found that i? 2 oo = 270±i 30 kpc, with the large spread 
coming primarily from the uncertainty in ot\. 
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Data Set 


V H 


Av* H 


r* H 

11 


Ar* H 

ll 


°8 


A<7 8 


t-test 


deBlok et. al. (53) 


119.0 


6.8 


3.62 


0.33 


0.613 


0.097 


1.36 


CF (348) 


179.1 


2.9 


7.43 


0.35 


0.84 


0.18 


0.43 


Mathewson et. al. (935) 


169.5 


1.9 


15.19 


0.42 


0.625 


0.089 


1.34 


Rubin et. al. (57) 


223.3 


7.6 


1.24 


0.14 


2.79 


0.82 


2.46 


Combined (1393) 


172.1 


1.6 


11.82 


0.30 


0.68 


0.11 


0.70 



TABLE I: The v* H (km/s), r* H (kps), and resultant a%, Aas, and t-test comparison with the WMAP 
value of as- The number of data points in each data set is listed in the parentheses. 

B. Estimating the Fractional Densities of Matter 



Because Pasymp is an asymptotic solution and has the same form irrespective the the detail 
shape of the galaxies, we can estimate fi asymp by averaging Pasym P ( r ) over a sphere of radius 
Xh/2 = r n , 



asymp 



(Pasymp)r J7 
Pc 



3Q> 



•57T 



1 + a A 

1 + 3a\ 



(l + «A) 2li+4 ] 



V(1+«a) / r \3 



(78) 

where p c is the critical density of the universe, and we have used Eq. (!70l) . Thus, the 
ratio O asymp /f2A depends only on the dimensionality and symmetry of spacetime, and the 
exponent a\. Numerically, we find f2 a sym P = 0.197±o.oi7- 

In performing this average we have implicitly assumed that there is only a single galaxy 
within the sphere, which is a gross under counting of the number of galaxies in the universe. 
Note, however, that p a sym P ( r ) is an asymptotic solution, and p}j(r) is a perturbation off 
Pasymp(^) that dies off rapidly with distance. While additional galaxies within the sphere 
may change the detail form of p aS ym P , these changes are expected to be equally short ranged; 
we thus expect Eq. ( 1781) to be an adequate estimate of fi asy mp- 

Such is not the case for f2 Dyn . Calculating f^Dyn directly by averaging p — p asymp would 
require knowing both the detail structure of galaxies, and the distribution of galaxies within 



S-0.025 
-0.026 



the sphere. Instead, we note that fl m = fi asy mp + ^Dyn, and using the value Q m = 0.239; 
from WMAP, find that fioyn = O.CMI^qq^. Thus, only a small fraction of the matter in the 
universe can be seen through their dynamics. 
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VI. CONCLUDING REMARKS 



Given how sensitive of our expression for <jg is dependent on v* H , r* H , and oca, that our 
predicted values of erg is within experimental error of its measured value is surprising. This 
is especially true as the data used in calculating <7g was taken by four different groups over a 
period of 25 years, and for purposes that have no connection whatsoever with our analysis. 
Even in the absence of a direct experimental search for a\, this provides a compelling argu- 
ment for the validity of our extension of the GEOM, and its impact on structure formation. 
This agreement also supports our free energy conjecture; the calculated of would be very 
different if j3 = 3, say, were used in calculating er| instead of (3 = 2. 

Direct detection and measurement of ct\ through terrestrial experiments may be possible 
in the near future. As mentioned in Sections II. D and III.B, at a value of 1.51 the exponent 
«a is likely small enough that the effects of the additional terms in the extended GEOM 
may soon be detectable. 

Interestingly, Vt m - Vt B = 0.196+o;q26 

is nearly equal to fi asymp in value. Correspondingly, 
Qb [i| is nearly equal to fi Dyn . It would be tempting to identify fi asymp with the fractional 
density of nonbaryonic (dark) matter in the universe, especially since matter in p asym p( r ) 
does not participate in the particle dynamics, and is not "visible" to measurements that 
inferred mass through particle motion under gravity. That ^Dyn would then be identified 
with Qb is consistent with the observation that most of the mass that has been inferred 
through gravitational dynamics indeed consists of baryons. We did not differentiate between 
normal and dark matter in our theory, however. Without a specific mechanism for funneling 
nonbaryonic matter into p asymp and baryonic matter into p — p asym p 5 we cannot at this point 
rule out the possibility that Q m — Qb = ^asymp and Qb ~ ^D yn is a numerical accident. 

APPENDIX A: FITS TO DATA 



In 



351 ]. measurements of the rotational velocity as a function of radius for 60 Sa, Sb and 



Sc spiral galaxies are given, allowing a fit of this data to v p tso (r). However, instead of fitting 



to t>P~ lso (r) directly as is done in 15], it is more convenient to fit the data to (t> p_lso (r)) 2 . 
As we are interested in the asymptotic velocity Vh instead of the density parameter for the 
pseudoisothermal profile, ours is a two parameter, (vh, Rc), least-squares fit to (t> p_lso (r)) 2 = 
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Vffc(r), where 

c(r) = 1 — arctan ( — — ) . (Al) 



It uses the variance 

N 

r 2 

(<uP-lSO)2 — J\f _ 2 



"~ — = ^ X>r°) 2 - ^c(r„)] 2 , (A2) 

n=l 

where {(v p _1so , r n )} is the set of velocity verses radius measurements for a galaxy with a total 
number of data points N. Minimizing with respect to v \ gives 



2 (« is °) 2 c(r n )) 
H «rnf) 



(A3) 



where (• • • ) denotes an average over the data points. Minimization with respect to Rq gives 
the implicit equation 

= <«-)M r „)) (|^J) - <*■.)•> (™) . (A4) 

Instead of solving Eq. (IA4I) directly, we substitute Eq. fl A3j) into Eq. (IA2I) . and find iteratively 
the i?c that minimizes cr 2 j . The value for v 2 H is then found through Eq. ( 1A3I) . 

Standard errors for vfj and Rc can be found directly. Taking the implicit derivative of 
Eq. (JAU), 

3R C _* c f<c W >? c(rj) /c^W\\ (A5) 



a(wf iso ) 2 A A \ rf + i? 2 , v v \ r 2 + i? 2 



where 



|fS>(f^>-«™-)KMw>- (A6) 

The standard error <j Rc in i?c is then 



For the standard error in vfj, we use Eq. ( 1A3I) and find 

gg| _ c(rQ ^ ^ 1 / <- is % 2 \ _ ( 2 ) / c(r n )r 2 \^ t| 



d(vf s y A(c(r n ) 2 ) V ^(c(r n ) 2 ) X^ + r 2 / Vn) 2 \i? 2 + r 2 /,/ i? c d(v^ 

(A8) 



to 
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resulting in a standard error in vh of 

1/2 



1 

2 



a (V iso )^ + A + 1 / vfr* \ 2 / c{r n )rl \ \ 2 a 2 Rc 



Nvl(c{r n y) V vl(c{r n y)\Rl + rl/ (c(r n )*) \B* + r*/ J R? c 

"(A9) 

n 

Our fits of the Rubin et. al. data are tabulated in Table HT1 The base data from [35[ was 
based on a Hubble constant of 50 km/s/Mpc, and the results given in Table HT1 are for this 

n 

value of the Hubble constant. Of the 60 galaxies from [351] , NGC 6314 and IC 724 could not 
be fitted to a nonzero Rc, while the fit for NGC 2608 resulted in a Rc that is less than 0.01. 

APPENDIX B: DATA SETS 

For the de Blok et. al. data set, the galaxy F568-3 was analysed twice; we use analysis of 



F568-3 given by the authors in |l5lj. In de Blok and Bosma 
and U5750, also appeared in [15J; we used the values from 



16J. two of the galaxies, F563-1 



15l | for these galaxies in our 



analysis. The radius for DD0185 from [15[] was not determined, and we could not include 
this data point in our analysis. Thus, out of 56 possible galaxies, 53 were used. 

For the Rubin et. al. data set, we could not find a nonzero radius for two galaxies, and 
one galaxy had a radius less than 0.01 kpc. As this radius was smaller than the resolution of 
their observations, this data point was not included. A total of 57 galaxies were thus used 
from [3^]. 

For the CF and Mathewson et. al. data sets, the vast majority of the data were based 
on single observations. We therefore had greater leeway in cleaning up this data, but even 



here we were circumspect. First, 75 galaxies in [36j were observed multiple times. Of these, 
the galaxies UGC 7234 and UGC 10096 had listed an asymptotic velocity for one of the 
observations that was opposite from the measured asymptotic velocity for the others. We 
assumed that this was a typographical error, and the sign of the rotational velocity for the 
anomalous velocity is reversed. Second, five galaxies in the CF and Mathewson et. al. data 
sets had a r# = 0, one galaxy had a radius core that was 11-sigma out, and three galaxies 
had a Vh that exceeded 8,000 km/s. These are likely indications that the data was not 
sufficiently accurate to allow for a fit of the velocity curve, and they were removed. Finally, 
given that there are only 1393 galaxies combined in the data sets, if a galaxy had a vh 
or a th that was six-sigma or more out from the mean, they were removed. In the end, 
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0.42 
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TABLE II: 



Fitted values of Rc (kpc) and vjj ( 
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(km/s), and their errors for Rubin et. al. data. 



Data Set 


Data Removed 


Reason 


-Data Set 


Data Removed 


Reason 


de Blok et. al. 


DD0185 


rn = oo 




ESO 243-G34 


vh is 5 a out 


Rubin et. al. 


NGC 6314 


r H = 




ESO 317-G41 


r H =0 




IC 724 


fH = 




ESO 358-G9 


vh is 6 a out 




NGC 2608 


r H < 0.01 




ESO 435-G25 


v h is 5 a out 


CF 


UGC 6534 


vh is 35 a out 




ESO 467-G12 


r H =0 




UGC 12543 


vjj is 11 a out 




ESO 554-G28 


u# is 6 <t out 


Mathewson et. al. 


ESO 140-G28 


v H > 8,000 km/s 




ESO 60-G24 


u/f is 10 a out 




ESO 481-G30 


v H > 24, 000 km/s 




ESO 359-G6 


rff is 11 a out 




ESO 443-G42 


v H > 94, 000 km/s 




ESO 481-G21 


r# is 6 o out 




ESO 108-G19 


r# = 




UGCA 394 


rjf is 7 a out 




ESO 141-G34 


vh = 




ESO 298-G15 


th is 7 a out 




ESO 21-G5 


r# is 6 <7 out 




ESO 545-G3 


rfl- is 7 cr out 




ESO 548-G21 


r# is 7 cr out 




ESO 404-G18 


r/f is 9 a out 




NGC 7591 


r H = 









TABLE III: Listed are the galaxies removed from the data sets used in our analysis along with the 
reason for their removal. 



348 galaxies were used in the CF data set, and 935 galaxies were used in the Mathewson 
et. al. data set. A summary of the data points not used in our analysis is given in Table HTT1 
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